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SYSTEM AND METHOD FOR ESTIMATING. SYNTHESIZING AND 



MATCHING NOISE IN DIGITAL IMAGES AND IMAGE SEQUENCES 

FIELD OF THE INVENTION 

5 The invention relates generally to the field of digital image 

processing and in particular to digital synthesis and matching of noise. 

BACKGROUND OF THE INVENTION 

It is generally known that compositing computer graphics and film 

1 0 imagery requires that the computer graphics elements be matched in terms of 
color and texture to the underlying background elements. Typically these 
background elements have been captured on film. It is also known that 
intercutting footage from different sources such as two kinds of film, or film and 
digital video, requires that the different elements match to create a seamless look 

1 5 for the entire image sequence. With more footage from mixed sources and 
computer graphics, there is a clear need for a robust grain simulation model. 
Moreover, there is a need for a fully automatic method of estimating grain 
parameters in image sequences to reduce the need for user interaction, which can 
be quite time consuming. 

20 Grain in digitized film images is well characterized by zero-mean 

Gaussian distributed noise for the three color channels. Zero-mean Gaussian 
models are fully characterized by their second-order statistics (the auto-covariance 
function). Modelling the interactions between color channels requires analysis of 
the spectral correlations. These spectral and spatial correlations can then be used 

25 to shape a white noise signal to create synthetic grain. In this way, synthesis of a 
simulated grain pattern is essentially a texture synthesis problem. But correct 
shaping of the white noise is required to generate a visually indistinguishable 
"match" for a given film grain. 

Texture analysis/synthesis methods are well known in the 

30 literature. Bergen and Heeger ("Pyramid-Based Texture Analysis/Synthesis" 

SIGGRAPH, 1995, pp. 229-238) describe a method for analyzing a texture pattern 
and then replicating it. This method suffers from improper color correlation 
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techniques generating significant bias in different color channels. The problem 
lies with the simple fact that once white noise fields have been convolved, the 
autocorrelation function (ACF) is non-zero for each color channel. Adding these 
together using a color transformation generates the correct spectral correlations 
5 but biases the spatial correlations. Biased spatial correlations will result in 
synthesized textures that look different from the desired texture. 

Chellappa and Kashyap considered autoregressive models for 
texture synthesis ('Texture synthesis using 2-d noncausal autoregressive models," 
IEEE Trans, on Acoustics, Speech, and Signal Processing, vol. 33, pp. 194-203, 

10 February 1985) but the optimization methods required in that work are inefficient 
(i.e., long processing times for both analysis and synthesis). 

The newer area of non-parametric synthesis requires templates of 
grain to perform the synthesis and is also known to be more computationally 
intensive than other methods. For example, Efros and T. Leung, ("Texture 

15 synthesis by non-parametric sampling," in ICCV99, pp. 1033-1038, 1999) 
considered how to sample pattern templates to generate new versions of the 
pattern. This is time consuming and cannot produce a wide enough variety of 
patterns when considering signal dependent noise. 

With regard to noise estimation, most methods use single frames. 

20 Consider, for example, Olsen ("Estimation of Noise in Images: An Evaluation," 
GMIP(55), No. 4, July 1993, pp. 319-323) and Immerkaer, ("Fast Noise Variance- 
Estimation", Computer Vision Image Under standing(64), No. 2, September 1996, 
pp. 300-302). These two methods estimate noise in generally flat regions of a 
single image. Motion compensated noise reduction methods are also known in the 

25 literature (see for example Dubois, E., Sabri, S., "Noise Reduction in Image 
Sequences Using Motion-Compensated Temporal Filtering," IEEE Trans. 
Communications (32), 1984, pp. 826-831). Although these methods are used for 
image sequences, they ignore the concept of iteratively estimating both motion 
and noise parameters. 

30 There are also relevant patents in the prior art. Published U.S. 

Patent Application 2002/0034337 (entitled "System for Manipulating Noise in 
Digital images") discloses a method for shaping white noise to create synthetic 
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noise with specific characteristics. The most obvious problem with this approach 
is the assumption that the color channels can simply be decorrelated and the noise 
synthesized independently. The subband decomposition procedure disclosed 
therein follows closely the treatment by Bergen and Heeger mentioned above. 
5 Indeed, it suffers from exactly the same improper color correlation techniques, 
thereby generating significant bias in different color channels. 

In commonly-assigned U.S. Patent No. 5,641,596 (which issued 
June 24, 1997 and is entitled "Adjusting Film grain properties in digital images" 
and), D. Cok and R. Gray disclose a method for synthesizing film grain using an 

10 autoregressive model that yields very good results. To achieve a better match to 
older or unknown film stocks, a larger spatial support is necessary. The model in 
that work had a limited spatial support and as scanning resolutions increase, 
spatial support needs to grow accordingly. 

The prior art disclosed above is deficient in its ability to efficiently 

1 5 generate a visually indistinguishable "match" for a given film grain. What is 

needed is a method of synthesizing noise that visually approximates a given film 
grain and is also close in a numerical sense to the statistics in a given random 
field. In particular, it would be desirable to perform grain matching when given 
the statistics for two different images. Finally, it would be desirable to have a 

20 robust method for automatically determining those statistics for image sequences. 

SUMMARY OF THE INVENTION 

The present invention is directed to overcoming one or more of the 
problems set forth above. Briefly summarized, according to one aspect of the 

25 present invention, the invention resides in a method and system for synthesizing 
noise in a digital image comprised of a plurality of color channels, where the noise 
is used to represent a synthetic texture that visually approximates a predetermined 
appearance. More specifically, the method comprises the steps of (a) estimating 
noise statistics based on a set of spatial autocorrelation coefficients and a set of 

30 spectral correlation coefficients that correspond to the color channels, where at 
least one of the spatial autocorrelation coefficients is weighted by at least one of 
the spectral correlation coefficients, thereby providing a weighted set of 
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autocorrelation coefficients, and where filter coefficients are obtained from the 
weighted set of autocorrelation coefficients; and (b) synthesizing a synthetic 
texture by using the filter coefficients to filter a random noise field in each of the 
color channels, thereby producing output noise fields in each of the color channels 
5 that replicate the synthetic texture when combined into the digital image. 

Accordingly, we disclose a method for synthesizing a texture when 
given second order noise statistics and spectral correlations. This data may be 
obtained from scanned film samples, from any digital image sequence, or may be 
interactively user controlled. We also disclose a method for matching noise 

10 between two different images using second order noise statistics and spectral 

correlations. Finally, we disclose a method for automatically computing required 
noise statistics in image sequences with or without test target data thus facilitating 
fully automatic robust noise matching. 

The principal advantage of the present invention is that it improves 

1 5 upon prior art by producing very high quality textures that resemble grain from 
film images. 

These and other aspects, objects, features and advantages of the 
present invention will be more clearly understood and appreciated from a review 
of the following detailed description of the preferred embodiments and appended 
20 claims, and by reference to the accompanying drawings. 



BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram illustrating the overall workflow for the 
method and system according to the invention. 
25 FIG. 2 is a block diagram illustrating a noise synthesis technique 

by which signal dependent noise is generated. 

FIG. 3 is a diagram illustrating a technique of filter coefficient 
generation by which target spatial autocorrelation parameters are modified by the 
spectral correlation coefficients to produce modified autocorrelation functions 
30 which are subsequently used to solve for the filter coefficients. 
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FIG. 4 is a block diagram illustrating a noise estimation technique 
by which signal dependent noise is estimated directly from image data, possibly 
utilizing an iterative step. 

FIG. 5 is a block diagram illustrating a noise estimation technique 
5 by which signal dependent noise is estimated directly from image data, possibly 
utilizing an iterative step, with the additional step of estimating and using spatial 
masks to improve the robustness of the noise estimation. 

FIG. 6 is a block diagram of an automatic noise-matching 
algorithm illustrating a technique by which the noise in a source sequence is 
10 matched to the noise in a target sequence. 

FIG. 7 is a block diagram of an automatic noise matching 
algorithm employing modification of multiple image sequences, which 
specifically illustrates a technique by which noise in one sequence is matched to 
the noise in another sequence by modifying both sequences to achieve a target 
15 noise. 

DETAILED DESCRIPTION OF THE INVENTION 

Because digital processing systems employing noise estimation, 
synthesizing and matching are well known, the present description will be directed 

20 in particular to attributes forming part of, or cooperating more directly with a 

system and method in accordance with the present invention. System and method 
attributes not specifically shown or described herein may be selected from those 
known in the art. In the following description, a preferred embodiment of the 
present invention would ordinarily be implemented as a software program, 

25 although those skilled in the art will readily recognize that the equivalent of such 
software may also be constructed in hardware. Given the system and method as 
described according to the invention in the following materials, software not 
specifically shown, suggested or described herein that is useful for 
implementation of the invention is conventional and within the ordinary skill in 

30 such arts. If the invention is implemented as a computer program, the program 
may be stored in conventional computer readable storage medium, which may 
comprise, for example; magnetic storage media such as a magnetic disk (such as a 
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floppy disk or a hard drive) or magnetic tape; optical storage media such as an 
optical disc, optical tape, or machine readable bar code; solid state electronic 
storage devices such as random access memory (RAM), or read only memory 
(ROM); or any other physical device or medium employed to store a computer 
5 program. 

The present invention improves upon prior art by producing very 
high quality textures that resemble grain from film images. Moreover, a technique 
for estimating the statistics from an arbitrary sequence is shown and a technique 
for matching the noise levels from two different sequences is given. 

10 Figure 1 shows the overall system and workflow for the invention 

as disclosed herein either for estimating the statistics from an arbitrary source or 
sequence (e.g., image sequence 1) or for matching the noise levels from two (or 
more) different sequences (e.g., image sequences 1 . . .K). In Figure 1, an image 
sequence 1 is shown having its noise reduced by a motion compensated filter 10 

1 5 followed by a noise estimation algorithm 20. That process may be performed 
once or may be repeated iteratively (30). At this point the image sequence on a 
line 32 has had its noise reduced and a good estimate on a line 34 of both the 
original and remaining noise is available. The image sequence may then have 
noise added back to it in an adder stage 40 in order, e.g., to visually approximate a 

20 desired appearance, such as the appearance of a given film grain. 

The characteristics of the noise that get added back may be 
determined in several ways from an arbitrary source or from one or more image 
sequences. First, the noise may be chosen from predetermined media "looks" or 
appearances 50 which will make the current sequence look, e.g., in terms of a 

25 given film grain, as though it originated on another media. This may be done by 
analyzing a given media or by having an artist design a particular look and then 
storing the parameters for that look. In this way, infinitely many media "looks" 
are possible. Alternatively, the noise may be manually controlled by a user in an 
interactive control system 60. Finally, the noise may be determined from the 

30 analysis of the noise estimate returned on line 34 from image sequence 1 or from 
an analysis (in an image sequence analysis stage 70) of one or more other digital 
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image sequences (e.g., the image sequence K). Any of these controls (50, 60, 70) 
can be used to set the "target" look. 

The noise matching block 80 performs the filter coefficient 
generation and the noise synthesis that determines the characteristics of the noise 
5 that needs to be added to the sequence or sequences in the adders 40, thereby 
forming an output image sequence 1 (and/or an output image sequence K) and/or 
a composite image sequence output 90. (In Figure 1, the additional "target" image 
sequence is shown as an image sequence K, which may be for example a sequence 
shot against a green screen that will be composited with the current image 
1 0 sequence 1 . The stage 70 for the image sequence K includes similar components 
10, 20, 30 and 40 for estimating and adding noise to its output image sequence K 
and/or to the composite image sequence output 90.) 

Signal dependent noise that very closely resembles film grain can 
be generated using the following model for generating noise: 

r(x,y) = a r r(x,y) 
g(x, y) - a g g{x, y) + f3 g r(x, y) 
15 b(x,y) = a b b(x,y) + (3 b g(x,y) Eq (1) 

This model is illustrated in Fig. 2 for a plurality of random noise fields in a 
plurality of color channels. In that figure, a random noise field 100 for the red 
channel is multiplied in a multiplier 110 by a scale factor and filtered in a filter 
120 to produce a scaled output noise field a r 7 (140). Similarly, a random noise 
20 field 102 for the green channel is multiplied in a multiplier 112 by a scale factor 
(Xg and filtered in a filter 122 to produce an intermediate noise field a g g . This is 

added (adder 131) to a portion f3 g of the red noise field 140, produced by a scaling 
multiplier 130, to generate a green noise field 142. Finally, a random noise field 
104 for the blue channel is multiplied in a multiplier 114 by a scale factor at and 

25 filtered in a filter 124 to produce an intermediate noise field a b b . This is added 

(adder 133) to a portion fib of the green noise field 142, produced by a scaling 
multiplier 132, to generate a blue noise field 144. The random noise fields 100, 
102, 104 are typically Gaussian random noise fields whose individual elements 
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are scaled using a noise table describing the variance at a given code value, i.e., at 
one or more signal levels. This noise table maps code values to variances or 
standard deviations, where the variances or standard deviations may be a function 
of the image data. Film grain is signal dependent so the underlying image content 
5 is needed to provide appropriate scaling factors (in this work we assume this 

scaling is done prior to scaling by the scale factors a r , a g , and a b .). It should be 

obvious to those skilled in the art that various combinations of the different 
channels may be added to yield a correlated color signal. 

Making noise appear as a predetermined texture, such as film grain 

10 requires that the filtering blocks 120, 122 and 124 shown in Fig. 2 produce correct 
target autocorrelation values and variances. In one aspect of the invention, a 
predetermined texture is synthesized when given a selected set of spatial 
autocorrelation coefficients p R (TV), p G (N), p B (AT) and a selected set of spectral 
correlation coefficients p RG , p GB , p RB . This coefficient data may be obtained from 

1 5 conventional sources, such as a calculation from scanned film samples or from 
any digital image sequence, or from interactive user control. Based on this 
coefficient data, it can be shown that the following relations will yield the correct 
noise statistics in the output noise fields 140, 142 and 144: 
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Pr(N) = p r (N) 

p s (N) = 



Pg (N) - P% P r(N) 



rg 



Pb(N) - p 2 gb p 9 (N) 



Eqs. (3) 

where the notation x(N) indicates a correlation at lag in one dimension. The 
equations are also valid when the lag is denoted x(m,n) where m and n denote a 
lag in the horizontal and vertical direction. This process of reweighting the 
5 autocorrelation coefficients is illustrated in Fig. 3 for a multi-channel source of 
noise data. Figure 3 illustrates how the input autocorrelation functions and 
spectral correlations 200 are used in the reweighting stage 210 according to Eqs. 
(3) to produce modified autocorrelations 212. Eqs. (3) describes the target 
autocorrelation values that the filtered noise fields should have so the final noise 

1 0 fields appear as specified when added together. Eqs. (2) describes how to modify 
the variance (the sigma terms) of the noise fields prior to filtering as well as how 
to weight the fields when creating the final synthesized noise (as shown in Figure 
2). With the target autocorrelation values (the tilde values) one can build a filter 
by invoking an optimization algorithm 220 for solving an optimization problem. 

1 5 The optimization problem solves for the filter coefficients when given the 

autocorrelation values. The formulation of the optimization problem depends 
upon the form of the filter. Filter coefficients for autoregressive models can be 
solved using the Yule- Walker equations; however, we use a modified linear 
system to minimize errors in correlation terms at higher lags. As an illustration, 

20 for terms up to lag 4 we can solve the following equation in a least-squares sense: 
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10 



- 10- 

The above steps are repeated for each color channel or band and for the horizontal 
and vertical image dimensions to build a separable filter. For autoregressive 
models the output is then obtained: 



ro = n ro + J2*=i Qih fo = r 0 + Y$Li fan 
90 = n go + Y%=i "4>i9i 90 = 90 + E?=i ^# 
bo = n bo + J2i=i Xi*>i b 0 = b 0 + ELi Xih 



Eqs. (5) 



If a moving average filter is desired, the filter coefficients may be determined by 
solving a nonlinear equation. An example for a one-dimensional non-linear filter 
with five coefficients is as follows: 

/ / m 0o$i+ 0102 + $263 + 9ah \ 2 , 

0 l A%A Wl V )fA) ~ el + 0l+ol + el + el ) + 

( ,„v 0002+0103+0204 ^ 2 , 

/ 0Q03 + 010* \ 2 , 



Eq. (6) 



Similarly, one may build a non-separable autoregressive or moving average model 
using the two dimensional autocorrelation values p R (m, ri), p c (m, n), p B (m, n) . 
These solutions are solved in a similar way to the methods disclosed above. 
Whatever the exact solution, the filter coefficients are clearly tied to the modified 

1 5 autocorrelation values. These filter coefficients are provided as the output 230 of 
the filter coefficient generation process, as shown in Fig. 3, and are then used as 
the filter coefficients for the filters 120, 122 and 124 in Figure 2. 

Sometimes the noise statistics of the source and target are known 
via a calibration procedure. More often than not, they are unknown and need to 

20 be estimated from data. Also, it should be clear that arbitrary noise matching is 
not always possible unless the source and target material are modified. Consider 
simply that one cannot add an uncorrected random signal to the source material 
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and end up with a lower noise position. If the target noise position is lower than 
the source position then some form of filtering is required. But noise filtering 
usually requires an estimate of signal noise which may be unknown. To 
accomplish noise estimation and filtering, Figure 4 shows a means of iteratively 
5 performing signal dependent motion compensated noise filtering and noise 

estimation. In Figure 4, a sequence of source images 300 is filtered by a filtering 
stage 310 and these filtered images are compared to the original images by a 
noise estimation algorithm 320 to produce an estimate of the noise 330. This 
process may be iterated as shown (332) where the new noise estimate is used in 

10 the motion estimation and filtering stage 310 to produce a better noise estimate. 

A preferred motion estimation technique for use with the present 
invention is disclosed in Black, M.J. and Anandan, P., 'The Robust Estimation of 
Multiple Motions: Parametric and Piecewise-Smooth Flow-Fields," Computer 
Vision and Image Understanding, January 1996, which is incorporated herein by 

1 5 reference. For the purposes of the present discussion a brief overview of the 

algorithm is given. Black et ah proposed the use of robust functions of the optical 
flow and smoothness terms in the general optical flow formulation. Their total 
energy function has the form: 

20 E t (t/, v) = X f p f (ul x + vl y + /, , o f ) + £ X s p s (u -u k ,o s ) + A s p s (v-v k a s ). Eq. (7) 

k 

The derivatives of the total energy function are: 



Be 

—^ = A f I x <p f (uI x +vI y +I,,a f ) + ^A s (p s (u-u k ,a s ) = 0 
-^7 = A/ 1 , <Pf « + vI y +/,,07) + £^(v-v A ,a-J = O 
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In these equations p is the robust error norm and q> is the derivative of p . For 
example, in Black et al. one robust function used is the Lorentzian 



p(x,a) = \og 



1+1 
2 

v 



X 



with ^ ) = _2fL_ 



5 The selection of the scaling term a in the above equations is 

varied according to the concept of graduated non-convexity, as described in the 
paper. However, in the present invention we propose tying the parameter to the 
signal dependent noise often found in scanned films. Doing this produces more 
robust motion estimates. 

10 Given a set of motion estimates, one may apply motion 

compensated noise filtering techniques, which are known in the literature to those 
skilled in the art. One means is disclosed in commonly-assigned U.S. Patent No. 
5,600,731 (which issued February 4, 1997 to M. Sezan, M. Ozkan, and S. Fogel, 
and is entitled "Method for Temporally Adaptive Filtering of Frames of a Noisy 

1 5 Image Sequence using Motion Estimation"), which is incorporated herein by 
reference. The filtered image is subtracted from the original image to produce a 
noise map and estimates of signal dependent variance and signal dependent 
autocorrelation parameters are computed using well-known statistical methods. 
These parameters may then be fed back into the motion estimator, as shown by an 

20 iteration operation 332 in Fig. 4, and used to re-compute the motion and filtered 
images. Such iteration produces more robust results. 

Occasionally the filtering produces artifacts, and the difference 
between the filtered image and the noisy image contains more than a 
representation of the noise but also signal content. To avoid biasing the statistics 

25 of the noise estimate, we avoid regions with large errors by using information 

from the filtering algorithm and/or information from the image itself to generate a 
spatial mask that may be used to modify the statistical computation of noise 
estimates. This procedure is shown in Fig. 5. As before, the entire procedure may 
be iterated (432) to improve robustness. In Figure 5, a sequence of source images 

30 400 is filtered in a motion compensated filtering stage 410 and these filtered 
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images are compared to the original images in a spatial mask generation stage 
420. A spatial mask is generated in the stage 420 to determine where good 
regions are in the image to estimate noise. The filtered images, original images, 
and masks are used by a noise estimation algorithm 430 to produce an estimate of 
5 the noise 440. This process may be iterated (432) as shown where the new noise 
estimate is used in the motion estimation and filtering stage 410, which produces 
more reliable masks and hence better noise estimate. 

The filtered images from a motion compensated noise filter can 
have noise added to them to produce matched noise where the matching is done 

10 relative to some other sequence of image data. This matching process is shown in 
Fig. 6. In Figure 6, target images 500 and source images 510 have their noise 
estimated by noise estimation algorithms 520 and 530. These noise estimates are 
used in a filter coefficient generation stage 540 to generate filter coefficients. 
Grain matching requires slightly different computations because the intrinsic noise 

1 5 present in the source and target image sequences needs to be handled. Eqs. (9) 
show the model where the subscript "T" denotes "target", "S" denotes "source", 
and "M" denotes "matching" 



Rt = Rs+PrRm 

Gt = Gs + PgGm + jgRm 

B t ^ Bs + PbBm+^bGm 

Eqs. (9) 

20 This model differs from our previous model since the source material has an 
intrinsic noise level. As before, we solve the weighting parameters: 
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Eqs. (10) 



Then the desired autocorrelation coefficients of the matching signal are: 
r (m,n) = — ~o Z2 



PG„(™,n) - 2 * 2 „^2 a 2 

a G T a Gs 7Q a R M 



PB M 



a B T ~ °B S ~ TB a G 



Eqs. 1 1 



5 The filter design process follows the exposition given earlier where either a 

separable or nonseparable filter may be used. With these filter coefficients, noise 
is synthesized in a noise synthesis stage 550 and added in an adder 552 to the 
source images resulting in "matched" source images 560 in which the noise 
position matches that of the target images. 

10 This filtering procedure may also be applied to both sequences first 

to lower the noise position and then different noise may be designed to give both 
sequences a similar noise position. Each source sequence may have a different 
source noise level. By choosing a similar target noise level, the matched noise 
added to each sequence will be different but the final noise positions will be the 

1 5 same. This method is shown in Fig. 7. In Figure 7, image sequence 600 

(sequence 1) is filtered by a motion compensated noise filter 602 and then the 
noise is estimated by a noise estimation algorithm 604, with this process repeated 
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iteratively (605) as desired. Filter coefficients are created in a filter generation 
stage 606 and noise is synthesized in a noise synthesis stage 608 and added to the 
filtered images in a summation stage 609 to produce a matched imaged sequence 
610 (sequence 1). 

5 Similarly in Figure 7, image sequence 620 (sequence 2) is filtered 

in a motion compensated noise filter 622 and then the noise is estimated by a 
noise estimation algorithm 624, with this process repeated iteratively (625) as 
desired. Filter coefficients are created in a filter coefficient generation stage 626 
and noise is synthesized in a noise synthesis stage 628 and added to the filtered 

1 0 images to produce a matched imaged sequence 630 (sequence 2) in a summation 
stage 629. The filter coefficients are tied together because the noise in each 
sequence is coupled to create the same "noise look" in both sequences. This is the 
same as the noise matching shown in Figure 1. Although the source noise level is 
different in image sequences 1 and 2, by manipulating the filter coefficients in 

1 5 generation stages 606 and 626 to achieve a similar target noise level, the matched 
noise added to each sequence in summation stages 609 and 629 will be different 
but the final noise positions of the sequences 1 and 2 will be the same. 
Thereupon, a conventional compositing operation 80 (Figure 1) may be performed 
on the matched sequences. The filter coefficient generation may be done 

20 automatically or may include some degree of manual intervention to guide the 
final look for each image sequence. 

The invention has been described in detail with particular reference 
to certain preferred embodiments thereof, but it will be understood that variations 
and modifications can be effected within the spirit and scope of the invention. 
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PARTS LIST 

10 motion compensated filter 

20 noise estimation algorithm 

30 iteration 

32 line 

34 line 

40 adder stage 

50 predetermined media "looks" 

60 user control system 

70 noise matching block 

80 composite output 

1 00 random noise field for red channel 

102 random noise field for green channel 

1 04 random noise field for blue channel 

110 red scaling factor multiplier 

112 green scaling factor multiplier 

1 1 4 blue scaling factor multiplier 

120 filter for red channel 

122 filter for green channel 

124 filter for blue channel 

130 red/green inter-channel scaling multiplier 

131 adder 

132 green/blue inter-channel scaling multiplier 

133 adder 

- 140 output red noise 

142 output green noise 

1 44 output blue noise 

200 input spectral correlations and autocorrelation functions 

210 reweighting stage 
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image sequence 2 
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622 motion compensated noise filtering 

624 noise estimation algorithm 

625 iteration 

626 filter coefficient generation 

628 noise synthesis 

629 summation stage 

630 matched image sequence 2 



